THE SMALL DEVIATIONS OF MANY-DIMENSIONAL DIFFUSION 
PROCESSES AND RAREFACTION BY BOUNDARIES 



Vitalii A. Gasanenko 

Abstract. We lead the algorithm of expansion of sojourn probability of many-dimensional 
diffusion processes in small domain. The principal member of this expansion defines nor- 
malizing coefficient for special limit theorems. 



Introduction. 

Let be a random process with measurable phase space (X, T,(X)). Consider the 
measurable connected domain D € and small parameter e. The investigations of 

asymptotics of sojourn probability (small deviations) 

P(f(t)eeD, t£[0,T]) (1) 

is jointed with many practice and theoretical problems [1-4]. In the literature, it was 
researched both rough asymptotics of principal member of (l)(log from it) [5] and exact 
asymptotics of diffusion processes of (l)[6-8]. In the works [9,10] was proved of algo- 
rithms of expansions of exact asymptotics of small deviation for diffusion and piecewise 
deterministic random processes for one-dimensional case. 

The purpose this article is to present the algorithm of expansion of small deviation for 
many-dimensional diffusion processes and to define all constants of principal member. 

In Section 1 our main result is stated and proved. In section 2 we consider the limits 
theorems about numbers of unabsorbed diffusion particles by boundaries of small domain. 

I. The expansion. 

We shall investigate of asymptote of following probability 

P{e,x) = P{£{t) e eD, 0<i<T), e -> 0, 
where G R d is solution of the following stochastic differential equation 

d 

dt(t) = a(t,Z(t))dt + Y,b i (Z(t))dw i (t), £(0)=*eeD. (2) 
i=l 

where functions 
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bi(x), a(t, x):R d ^ R d and R+ x R d -> R d . 

are differentiable. 

Set a^s) =£&£(*)&£(*). 

k 

It is known that P(e, x) = u^(T, x). Here u^t, x) is solution of the following parabolic 
boundary problem at < t < T 

i,j=l % 3 i=l 1 

u(t,x)\ t =Q = 1; xeD e ; u(t,x) = x G dD e , < t < T. (3) 

where D e = eD. It is assumed that D is a connected bounded domain from R m ; the 
boundary dQ is the Lyapunov surface C^' A ^ and Oefl. We interest of the asymptotic 
expansion e — > of solution this problem Uq(£, x) at e — > 0. 

We define the differential operator A : ^ a ij (0) dx dx ■ Let ^ be a matrix with 

i<i,j<d 1 J 

the following property 

^2 a lj (0)z l z j > yu|z| 2 . 

l<i,j<d 

Here /i, there is a fixed positive number, and z = (z±, ■ ■ ■ , zj) is an arbitrary real 
vector. 

This operator acts in the following space 

H A = {u:ue L 2 (D) nAue L 2 (D) n u(dD) = 0} 

with inner product (u 1 v)a = (Au,v). Here (, ) is inner product in L 2 (Q). The opera- 
tor A is a positive operator[ll]. It is known that the following eigenvalue problem 

Au = -Xu, u(dD) = 

has infinite set of real eigenvalues Ai — ► oo and 

< Ai < A 2 < • • • < A s < • • • . 

The corresponding cigenfunctions 

fill • • • i flni ; ' ' ' i fsl ; ■ • ■ j fsn s 7 

form the complete system of functions both in Ha and L 2 (Q) := {u : u G L 2 (Q) n 
u(dQ) = 0}. Here the number nk is equal to multiplicity of eigenvalue Afe. 

It is often convenient to present the system of eigenfunctions by one index: {/„(z)}. 
The corresponding system of eigenvalues {A„} will be with recurrences. We shall use it 
too. 

We introduce the spectral function 

e(x,y,X) = fi( x )fj(y)- 
\j<\ 

We shall need in the following theorem from the monograph [12]. 
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Theorem 1 ([12]. Th. 17. 5. 3). . There exists such constant C a that 



sup J\D« y e(x,y,\)\<C a \( n +W 2 

x,y£D v 

Here a is multi-index. 
Theorem 2. .If the surface 3D is Lyapunov surface and 



( dai(t,z). ,dbi(z) dai(T-t,z).\ 
s U p max I'l-flT-l'l 7ft 1 

T]xD,l<i,j<d (. oz 3 oz ] 01 ) 



< OO 



(t,z)e[0,T]xD,l<i,j<d 

then the following relation takes place at e — > 



P(e,ze) =exp|-Ai-2 +J n(t)dt\ ^ c lm f lm (z) (I + O(e)) , at z e D, 



where 



and ci m = / f lm (z)dz. 

D 

Proof. Make the change of variables and function 



= z*e, ^^exp je^a fe (T-t, 0)z fc j ' 



Now we obtain the following parabolic problem for function u\ 



duUt,z) 1 v-^ , .d 2 uUt,z) lv^/ _ . 1 . . , .\ duUt,z 



^ 2 e 2 .4-1 ' 0*0^ e^V 2 -jv-/-jv- ■'"') dz, | 

2 ( - -*»0) - ««Oi(T - t,0) Oj -(T - i,e*) - e ^ 2 /" — *») «i, « G I>;| 



ul(t,z)\ t =o = exp|e^a fe (T,0)^| ; z£fl; u e 1 (t,z) = z e 9-D, < t < T. 

(4) 

We will construct the asymptotic expansion of solution for this initial - boundary 
problem in the following form 

*4(t,z) = ^ fe (i,z)e fc . (5) 

fc>0 

Note that the famous expansion 
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cx P Lj2ak(T,0)z k \ = l + e^a fe (T,0)z fc + i ( e ^ a k (T, 0)z k J + ••• , 
I fe=i J fc=i ' \ fc=i / 



defines the initial conditions for v k , k > 0: 



d 1 / d \ 2 

v (0,z) = l, v 1 {0,z) = ^2a k (T,0)z k , v 2 {Q,z) = -\S^a k {T,Q)z k \ ■ ■ ■ . 

k=l \fe=l / 

Using the first fragment of Taylor series in zero point under conditions of theorem we 
can obtain the following representations 



a y - {ez)=a ij {0) + ea € ij (z), a, (T - t , ez) = a t (T - t , 0) + ea* (T - t, z) , 1 < i, j < d (6) 
where 



sup |af-(z)| < oo, sup \a\{T — t,z)\ < oo 

zeD,ee[oa]a<ij<d zeD,te[o,T],ee[oA]A<i<d 

Now, after substitution of (5), (6) to (4) we conclude that the Vo satisfies the problem 
v \dD = 0; v (0, z) = l, z e D. 

Here 

Kt)=^(^ l3 (0)a l (T-t,0)a 3 (T-t,0)-S lJ a l (T-t,0)a J (T-t,0)] . 

Further, let us denote by B e (t,z) the operator C 2 (D) -> C(D), for / G C 2 (D) it's 
defined as follows: 

B f (t,z)f = 

1 \^ f / x 9 2 f lr^/ , m , 1 , , /m „.\9/ 

& E ^w&^ + 7£r r "* ,B) "2 ff « {tt)flj(r - t,0) ^ + 
+e E (^«w°*( T - *>°K( T - *.°) - <m t - *>°) a K T - *. z ) - aat(T g^' 0) ^ 
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Now, formally the functions v k , k > 1 are denned by the following recurrence system 
problems 



Ov 1 I 



(8) 



w o |an = 0; w fe (0,z) = ^^a fc (T-i, 0)z fe ^ , 



We shall solve the problems of (7), (8) by method of separation of variables. According 
to this method the solutions are defined in the form 



v k (t,z) = ^q k . n (t)f n (z). 



(9) 



n>l 



For definition of principal number it suffices to construct of the vq. If we substitute 
(9) at k — to (7) then we obtain 

^2 | _ 9o,n(i) - ^r<7o,™(£) +M*)<?o,n(*) j fn(z) = 0. 

n>l ^ ' 

Set 

co,n — J £> fn(z)dz (coefficients of expansion of indicator of set D). The initial 
condition of vq has the following stating 



vo 



(0, z) = <?0,n(0)/„(z) = Co, n f n (z) = ^ X! c 0,lmflm(z), Z <E D. 



n>l 



n>l 



2>1 m=l 



By definition of system of functions {/ ra (z)}, now we have the system of ordinary 
differential equations 



qO,n{t) + ( -T - ) 90,n(i) = 0, go,n(0) = Co, n - 



From the latter one we have 



Set 



qo,n(t) = c ,„exp j-^T* + J ^(s)ds^ . 
A = sup \alj(z)\, L = ^2 (co.mi) 2 ■ 



e<l,z£D;i,j 



A\ = sup 

0<e<l,zeD,tE[0,T];i,j 



Z>l,l<m<n; 



di(T - t,ez) - -o-ij(ez)a,j(T - t,0) 



sup 

0<e<l, z£D,t£[0,T];i,j 



l -o\ 2 (z)a t (T - t,0) aj (T - t,0) - Si jai (T - t,Q)a){T -t,z)- da:> 



dt 



We have the following relations for eigenvalues A; 
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kil 2/d <Xi< k 2 l 2/d , max(fci, k 2 ) < oo 
Applying Cauchy-Bunyakovskii inequality, Theorem 1 and the latter one, we get 



E«k(*) 



d 2 v 



dzidzj 



^cxpf-A ; ie 2 + / /z(s)ds^ ^V^ijW 
i ^ Jo / m=1 i • 



d 2 / m ,M 



dzidzj 



dzidzj 



< 



< 



< Aod^exp (^-\it€- 2 + ii{ s )d^j (Jr (c . ml ) 2 ^ ^jr ^ 

< A dC 2 , 2J LoEexp + J ^{s)d)j \f +2 < cxp (~^r) #0- (10) 

Here ifo < oo. 

Reasoning similarly we convince ourselves that for other parts of B e (t 7 z)v n the fol- 
lowing estimations take place 



A{(t,z)vo\ < AidCi^Lp y^exp ( + f fj,(s)ds J A ; 2+1 < exp 

«>1 \ £ n / 



¥) (ID 



|A 2 (t,z)«o| < A 2 dCo, L J2 cx V + J Af < exp (-^) #0,2, (12) 

where max{_Ko,ij -K0.2} < 00. 

Now let us estimate the coefficients f3„(t) of expansion of B e (t, z)vo by system {f n }n>i- 
Applying (10)-(12) and Cauchy-Bunyakovskii inequality, we get 

= I J B*(t, Z )v (t,z)f n (z)dz\ < 
D 



< 



J(B c (t,z)vofdzj yf f 2 n {z)dz\ < 



< exp(-A 1 te- 2 ) ( K " ' Au ' : ■ ,/v,,. 



The latter one now gives 



I J P e n (s)ds\ < e 7£ (t), 




(13) 
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where 

sup 7 e (i) < oo. 

0<e<l,te[0,T] 

Finally, let us estimate the difference r £ (i, z) = (t, z) — Vo(t, z). By definition, r e (t, z) 
is solution of the following problem 



dr 



l = ^^a«W^y + B e {T-t t z)vo zeD; (14) 



r e {t,z)\ t=0 = cxp je^a fe (T,0)z fc j - 1; zefl; r e (i,z) = z G dD, <t <T. 

It is clear that r £ (0, z) we can present as er{(0, z), where r\ (0, z) is uniform bounded 
function of variables e G [0, 1] and z G D. So, the coefficients of expansion this function 
by system {/„(z)} have the following forms 

/ r £ (0, z)f n {z)dz = e^ n , where sup V (^) 2 = M < oo. (13) 

Now we have the solution of (14) in the following form 

r e (t,z) = eJ2^ n cxp{-X n te- 2 + / («)ds}/„(z) 



n>l 



Applying latter one ,(13), (15), Theorem 1 and Cauchy-Bunyakovskii inequality we get 
at t > 

\r e (t,z)e~ 1 \ < (j>£) 2 ] C ,o^exp{-A^e- 2 + f F n (s)ds}\l} < 

\n>l J n>l J 

< MCq,o exp{ — \ite~ 2 }K ,3, where K n _ 3 < oo. 
The proof of theorem is completed. 

Remark 1. According to the above system of problems for definition of the functions 
Vk,k > 1, we outline the construction of coefficients qk,n)(t) for the series (8): 

ikA 1 ) = + {jk + Mfe-i,n(*)) q k . n (t), 

9fc,n(0) = J v k (0,z)f n (z)dz = i J ^2 a m (T,0)zSj f n (z)dz 



Here H%_ hn {t) = J f n (z)B e (t, z)v k -i(t, z)dz. 

D 
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Remark 2. Theorem 2 is coordinated with results of works [6-8] where the principal 
member of small deviations in ball are investigated for more simple SDE. 

II. The rarefaction of set of diffusion processes by boundaries of small 
domains. 

The following problem was investigated in works[13,14]. Let a set identical diffusion 
random processes start at the initial time from the different points of domain D. These 
processes are diffusion processes with absorbtion on the boundary dD. We are interested 
in distribution of the number yet absorbed at the moment T. The initial number and 
initial position of diffusion processes are defined either a random Poisson measure[14] 
or deterministic measure [13]. The proved limits theorems described the situation when 
T — > oo and initial number of diffusion processes depended on T and it increased at the 
rise of T. The role of normalizing function played principal member of asymptote of 
solution of according parabolic problem at T — > oo. 

Henceforth we shall assume that considered diffusion processes satisfy of the SDE (2) 
with different initial points. 

Now we consider the situation when initial number of absorbing diffusion processes 
in small domain eD depends on e — > and it increase under the condition of decrease 
of e. It is not hard to show, that now normalizing function is the principal member of 
parabolic problem (3) at e — > 0. 

The proofs of stated below theorems repeat the proofs of according theorems from 
[13,14] almost word for word. 

We will denote by rj(e, T) the number of remaining processes in the region eD at the 
moment T. 

We will also assume that tr-additive measure v is given on the E„- algebra sets from 
D, v{D) < oo. All eigenfunctions : D — > R 1 are (E„,Ey) measurable. Here Ey is 
system of Borel sets from R 1 . Let => denote the weak convergence of random values or 
measures. 

At the beginning we assume that initial number and position of diffusion processes are 
defined by deterministic measure N(eB, e), B 6 D. Thus, N(eB, e) is equal to number of 
starting points in the set eB. 

Let us denote by v e (-) the measure 



Theorem 3. Under the assumptions of the Theorem 2 let the N(e-,e) satisfies the con- 
dition 




where B € E„ . 

By definition of measure v t { ), we have 




exp ( 



TAi 
"2?" 




otherwise. 
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Then rj(e, T) =>■ rj(T) if e — > where rj(T) has Poisson distribution function with 
parameter 



a(T) = cxp | J fj,(s)ds J J F{z)dv{z), 



where F(z) = J2 fu{z)c u , c u = J fu(z)dz 

»=1 Q 
and n(t) is the function from Theorem 2. 

Now we consider the case when the initial number and positions of processes are 
defined by the random Poisson measure ^(-, e) in eD: 

P { ^eA,e) = k) = ^^e-^\ 

where m(e •, e) is finitely additive positive measure on eD for fixed e. 
We assign 



5(e) = cxp (-T^r) • 



Theorem 4. Under the assumptions of the Theorem 2 we suppose that m(e-,e) holds 
the condition 

lim m(eB, e)g(e) = v(B), B e £„. 

€— >0 

Then r](e, T) =^> ?7(T) «/ e — > w/iere ?7(T) /ias i/ie Poisson distribution function with 
the parameter a(T) from Theorem 3. 
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